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STRESS  AND  VELOCITY  FIELDS  IN  SOIL  MECHANICS1 
by  R.  T.  Shield 
Brown  University 


Summary 

Lines  of  discontinuity  in  the  plastic  stress  field  of 
a cohesive  soil  in  a state  of  plane  strain  are  discussed*  The 
Jump  conditions  on  the  stresses  and  the  restrictions  on  the 
velocity  field  in  the  neighborhood  of  the  line  are  obtained* 

The  theory  is  applied  to  the  problem  of  uniform  pressure  on 
one  face  of  a wedge  (or  earth  dam) , and  to  the  problem  of  a 
loaded  trapezoid. 

Finally,  an  analytic  integration  of  the  plane  strain 
equations  is  carried  out.  Expressions  are  derived  for  the 
coordinates  and  curvatures  of  the  failure  lines  and  the  velocity 
components  at  any  point  of  a plastic  stress  field  in  terms  of 
the  boundary  values* 


The  results  presented  in  this  paper  wore  obtained  in  the  course 
of  research  sponsored  by  the  Office  of  Naval  Reuearch  'under 
Contract  N7onr-35801  with  Brown  University. 
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i.  intros  U<?n» 

The  solutions  of  many  two-dimensional  problems  in  soil 
mechanics  are  based  upon  the  assumption  that  the  soil  is  a plastic 
material  in  which  slip  or  yielding  occurs  when  the  stresses 
satisfy  the  Coulomb  formula  Cl] 

2 2 \1/2 

f = sin  9 + £f^dx"’ °y^  + x xyj  ~ c cos  9 = 0,  (1) 

where  c is  the  cohesion  and  9 is  the  angle  of  internal 
friction  of  the  soil.  Neglecting  the  weight  of  the  coil,  the 
stresses  also  satisfy  the  equations  of  equilibrium 


0o_  0T: 

Tx  * “fly*  * °> 


l 


ax 


(2) 


The  two  characteristic  lines  of  the  hyperbolic  system  of  equations 
(1),  (2)  are  inclined  at  an  angle  itA  + 9/2  to  the  direction 
of  the  algebraically  greater  principal  stress  at  any  point.  As 
in  [2]  the  characteristic  linos  will  be  called  the  first  and 
second  failure  lines,  with  the  convention  that  the  direction 
of  the  first  failure  line  at  a point  is  obtained  from  the 
direction  of  the  algebraically  greater  principal  stress  by  a 
clockwise  rotation  of  amount  *A  + 9/2.  The  inclination  of  the 
first  failure  lino  to  the  x-axis  will  be  denoted  by  0. 

Denoting  the  principal  stresses  by  d^  and  <*2^dl  < d2^» 
it  is  convenient  to  introduce  the  quantity 

(d2-V 

p " 2 sin  9 ^ ^3) 

^Numbers  in  square  brackets  refer  to  the  bibliography  at  the  end 
of  the  paper. 
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An  arbitrary  constant  which  has  the  dimensions  of  stress  is  some- 
times introduced  -into  the  denominator  of  the  expression  for  p 
in  order  to  make  p non-dimensional,  but  this  practice  will  not 
be  followed  here.  With  the  yield  condition  (1),  it  can  be  shown 
[ 3 ] that 


<JX  = -p  [l  + sin  <p  sin(2©  + 9)1  + c cot  9, 

0y  = -p[l  - sin  9 sin(2©  + 9)1  + c cot 
axy  ” P s*n  V cos(2©  +■  9), 

and  the  equations  of  equilibrium  (2)  can  be  replaced  by  the 
equations 


^ cot  9 log  p + © = const,  along  a first  failure  line,  "| 
cot  9 log  p - © = const,  along  a second  failure  line, 


Of) 


(5) 


These  equations  were  obtained  by  Massau  [4l  for  a cohesionless 
soil  and  by  Kotter  15  J for  a cohesive  soil, 

Drucker  and  Prager  [6 1 considered  a proper  generalization 
of  the  Coulomb  hypothesis  (1)  and  used  the  concept  of  plastic 
potential  [ 7l  to  obtain  a stress-strain  law,  assuming  that  the 
soil  is  a perfectly  plastic  body.  In  the  case  of  plane  strain, 
the  stress-strain  law  corresponding  to  the  yield  function  (1)  is 

Sf 

x ' r 

where  ex,  e^,  y^  aro  the  plastic  strain  rates  and  ^ is  a 
non-negative  factor  of  proportionality.  With  equations  (4-),  the 
relations  (6)  can  be  written 


\ df 
y 


Yxy  ~X'd‘ 


(6) 
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ex  B jfsin^  - sin(20  + 9)}  , 1 
e » ^{sin  * + sin(2y  t 9)}  , 

Yxy  = X cos(26  + 9). 

It  will  be  assumed  that  there  is  no  deformation  of  the  soil 
until  plastic  yielding  occurs  and  therefore 


(7) 


8v  _ 
ey  " ey» 


Hi  + fcv 
6 y Ox 


» 


where  u,v  are  the  components  of  velocity  along  the  x,y-axes. 
The  rate  of  dilation  is  found  to  be 


ex  + ty  « X sin  9^0,  (8) 

so  that  plastic  deformation  must  be  accompanied  by  an  increase 
in  volume  if  9 / 0. 

It  was  shown  in  [2]  that  the  relations  (7)  imply  that 
the  rate  of  extension  along  the  failure  lines  is  zero.  Also, 
the  characteristic  lines  of  the  velocities  coincide  with  the  char- 
acteristic lines  of  the  stresses  and  it  is  therefore  convenient 
to  refer  the  velocity  equations  to  the  characteristic  lines.  The 
orthogonal  projections  of  the  velocity  vector  on  the  directions 
of  the  first  and  second  failure  lines  passing  through  the  point 
are  denoted  by  v^  and  Vg.  The  velocity  projections  v^,  Vg 
are  related  to  the  cartesian  components  u,v  of  the  velocity 
by  the  equations 

v^ = u cos  0 + v sin  0,  V2  « -u  sin(0+9)  + v cos  (©49),  j 


u =[vic°s(0+9)-v2Sin  0}sec  9,  v efv^sinC©^  )+V2Cos  Gjsec  9,  J 


(9) 


I t 

I 
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The  condition  that  the  rate  of  extension  along  the  failure  lines 
is  zero  is  expressed  by  the  equations 

dv,  - (v-jtan  9+  V2  sec  <p)d0  = 0 along  a first  failure  line, 
dv^  + (v^sec  9+  v2  tan  9)d0  = 0 along  a second  failure  line,^ 

and  these  equations,  together  with  the  velocity  boundary  condi- 
tions and  the  condition  of  non-negative  dilatation,  determine 
the  velocity  field  when  the  failure  lines  are  known. 

Discontinuities  in  the  velocity  field  were  considered 
in  [2]  and  it  was  shown  that  a line  of  discontinuity  in  the 
velocity  field  must  be  a failure  line.  The  change  in  velocity 
across  the  line  must  be  inclined  at  an  angle  9 to  the  line  of 
discontinuity,  so  that  a discontinuity  in  the  tangential  velocity 
is  accompanied  by  a separation  or  a discontinuity  in  the  normal 
velocity.  In  terms  of  the  velocity  projections  v^  and  v>>, 
the  requirement  that  the  change  in  velocity  is  inclined  et  an 
angle  9 to  the  line  Implies  that  either  or  V2  is  con- 

tinuous across  the  line,  according  as  the  line  is  a second  or 
first  failure  line  respectively. 


In  the  theory  of  a Prandtl-Reuss  material,  for  which 
9 - 0,  the  possibility  of  lines  of  discontinuity  in  the  stress 
components  is  well  known,  and  discontinuous  stress  fields  have 
proved  useful  in  many  problems.  Discontinuous  stress  fields  are 
of  value  when  limit  analysis  is  used  to  obtain  lower  bounds  for 
the  collapse  values  of  the  surface  tractions  in  a body  or 
assemblage  of  bodies,  even  though  the  stress  fields  may  be 
without  obvious  physical  significance.  It  seems  worthwhile, 
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therefore,  to  consider  discontinuous  plastic  stress  fields  in  a 
soil,  not  only  from  the  viewpoint  of  obtaining  real  solutions 
but  also  with  regard  to  future  applications  of  limit  analysis 
to  soil  mechanics*  A discussion  of  stress  discontinuities  in  a 
cohesionless  heavy  soil  has  been  given  previously  by  Massau  [4]. 
The  following  discussion  applies  to  a cohesive  soil. 

In  Fig.  1,  the  line  -t-'t  represents  an  element  of 
a line  of  stress  discontinuity  separating  the  plastic  stress 
fields  3 and  Jj.  Subscripts  a and  will  be  used  to 
distinguish  the  values  which  a quantity  assumes  on  the  two 
sides  of  the  line.  The  normal  to  the  line  is  inclined  at  an 
angle  Q to  the  x-axis  and,  from  equations  (4),  the  normal 
and  tangential  tractions  on  such  a line  are 

N = -p  [1  + sin  9 sin(20-2Q+9)  ] + c cot  9, 

T = p sin  9 cos(2G-2Q*-9). 


(11) 


The  equilibrium  of  the  small  rectangular  element  shown  in  the 


<•<  requires  that  the  normal  and  tangential  components  of 


* A f,  4AA 


stress  N and  T are  continuous  across  the  line,  but  the 


interior  components  Na, 


N*  may  be  discontinuous.  From  equa- 
b 


tions  (11),  the  equilibrium  conditions  are 


Pa[l+sin9  sin(2©&-2Q+9)3  = pbCl+sin  9 sin( 2^-204-9) ] , 

* 

pa  cos(2«a-22+9)  = pb  cos(29b-2Q+9). 

The  elimination  of  p , p between  these  two  equations  gives* 

a 0 1 

e.fter  some  reduction, 


sin(©a+©b-2af9)  + sin  <p  cos  (©a-©b)  = 0, 


(13) 
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provided  that  sin(©a-0^)  / 0,  that  is,  provided  that  the  dis- 
continuity is  not  of  zero  strength# 

The  condition(13)  together  with  one  of  the  conditions 
(12)  correspond  to  the  Jump  conditions  established  by  Prager  [8] 
for  a Prandci-Reuss  material  (9  = 0).  For  zero  angle  of  friction, 
condition  (13)  becomes 

©a  6b  = n7U  (1*+) 

The  Mohr's  circles  for  the  regions  a and  £ are 
shown  in  Fig*  2,  where  p&  is  taken  greater  than  pb  for 
definiteness.  The  circles  touch  the  yield  locus 

|t|  = c - o tan  9, 

since  the  regions  are  at  the  point  of  yielding.  The  stress  point 
C has  the  coordinates  (N,  -T)  and  the  circles  intersect  at  this 
point.  The  poles  of  the  two  circles  are  obtained  by  drawing  a line 
through  C parallel  to  the  element  t - of  the  line  of  dis- 

continuity. The  points  Pft,  P^  where  this  line  meets  the  circles 
are  the  poles  of  the  circles.  If  and  A^,  3^  denote 

the  points  of  contact  of  the  circles  3 and  Jfc  with  the  yield 
locus,  then  the  lines  P^,  PaBft  and  PbAb,  P^B^  give  the 
directions  of  the  first  and  second  failure  lines  in  the  regions 
a and  £ respectively.  It  can  be  seen  from  the  diagram  that 
the  direction  of  the  line  of  discontinuity  lies  in  the  acute  angle 
formed  by  the  failure  lines  in  region  £ and  in  the  obtuse  angle 
formed  by  the  failure  lines  in  region 

It  can  easily  be  shown  that  a failure  line  cannot  be 
a line  of  discontinuity  in  the  stresses.  Since  a discontinuity 
in  the  velocity  field  can  occur  only  across  a failure  line,  it 
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follows  that  the  velocity  field  must  be  continuous  across  a line 
of  stress  discontinuity*  For  a straight  line  of  discontinuity, 
which  may  be  taken  to  be  parallel  to  the  x-axis,  the  continuity 
cf  the  velocity  component  u across  the  line  implies  that 
8u/8x  = ex  is  continuous  across  the  line.  This  argument  can 
be  extended  to  show  that  ex  is  continuous  across  a curved  line 
of  stress  discontinuity  at  any  point  P of  the  line,  where  the 
x-axis  is  taken  parallel  to  the  tangent  to  the  line  at  P. 

Considering  first  a Prandtl-Reuss  material,  equations  (7) 


give 


sin  26, 


since  9 is  zero  for  this  material.  It  follows  that 


\a  sin  20a  sin  20b,  (1?) 

where  X^,  ^ and  6&l  ©b  are  the  values  assumed  by  X and  0 
in  the  two  stress  fields  at  the  point  P.  The  Jump  condition  (1*+) 
on  0a,  ©b  requires  that 

©b  " -®a  - (16) 

Since  the  line  is  not  a failure  line,  sin  28a  t 0,  and  the 
substitution  of  (16)  in  (15)  gives 


Ka  = ~V 


This  equation  implies  that  X 
point 


and  X.  are  both  zero  at  the 


a — b 

P because  Xft  and  X^  are  non-negative  quantities.  The 
plastic  rate  of  strain  is  therefore  zero  at  points  on  the  line  of 
discontinuity,  and  the  line  must  be  considered  as  a filament  of 
non-plastic  material.  Elastic  strains  are  neglected  so  that  the 
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filament  must  be  taken  to  bo  inextensible  but  perfectly  flexible. 
This  conclusion  has  been  reached  previously  [9,10]  by  assuming  that 
the  line  of  discontinuity  is  the  limiting  case  of  a narrow  transi- 
tion region  in  which  the  stress  varies  in  a rapid  but  continuous 
manner. 

In  the  same  way,  a line  of  stress  discontinuity  in  a 
soil  must  be  considered  as  an  inextensible  filament  of  non-plastic 
material.  Apart  from  the  algebra,  the  discussion  is  the  samo  as 
the  above  discussion  for  a Prandtl-Reuss  material,  and  it  will  not 
be  given  here. 

3*  Wetas  wter  gnUatasttl  Rreaavygt 

The  problem  of  wedges  with  uniform  pressure  on  one  face 
will  be  considered  in  this  section  and  in  the  following  section. 

A continuous  stress  solution  for  obtuse  angled  wedges  has  been 
given  by  Prandtl  [ll]  and  the  pattern  of  the  failure  lines  is 
shown  In  Fig.  3a.  The  wedge  ABE  of  angle  Pq  (Pq  n/2 ) is 
loaded  by  a uniform  pressure  P along  AB,  producing  a plastic 
state  of  stress  in  the  region  ABEDC.  The  regions  ABC,  BDE  are 
regions  of  constant  stress  and  the  region  3CD  is  a zone  of  radial 
shear  of  angle  Pq-V2.  The  pressure  P required  to  produce 
plastic  flow  with  this  pattern  of  failure  lines  is  given  by 

P = c cot?  £exp[(2pQ  - r.)  tan  ? ]tan2(  ^ ) - l|  • (17) 

The  stress  field  of  Fig.  3a  is  not  applicable  to  acute 
angled  wedges  (p^  < n/2)  since  the  two  constant  stress  regions 
ABC,  BDE  would  overlap.  However,  the  introduction  of  a line  of 
stress  discontinuity  in  the  plastic  stress  region  enables  a stress 
solution  to  be  found.  Fig.  3b  shows  the  wedge  loaded  along  AB. 
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The  constant  stress  regions  ABC,  BCD  are  separated  by  a line  of 
stress  discontinuity  BC  which  is  inclined  to  AB  at  an  angle 
Y to  be  determined.  AC  and  CD  are  first  failure  line3.  The 
different  values  a quantity  may  assume  in  the  regions  ABC,  BCD 
will  be  distinguished  by  the  subscripts  a and  respectively. 
Choosing  the  x,y-axes  as  shown  in  the  figure,  the  inclination  0 
of  the  first  failure  lines  assumes  the  values 

©a  e °>  ®b  = -(n/2  ‘ P0)# 

The  substitution  of  these  values  into  the  jump  condition  (13), 
in  which  Q is  put  equal  to  y-  (*A  - 9/2),  gives  the  relation 

sin  (P0  - 2y)  4 sin  9 sin  = 0. 

This  equation  determines  the  angle  Y and  the  relevant  root  of 
the  equation  is  found  to  be 

y = Pq/2  ♦ ji/2,  (18) 

where  n is  given  by 

sin  p = sin  9 sin  0 £ p.  £ */2.  (19) 

The  value  of  p in  the  region  BCD  is  determined  by  the  condition 
of  zero  traction  on  BD  and  it  has  the  value 

c cot  9 

pi,  = 9 

b (1-sin  9) 

The  second  equation  of  equations  (12)  then  gives 

- c cot  9 sinffip  •»  p.) 
a (1-sin  9)  sin(p0+ii) 


I 
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With  equations  (11),  the  normal  pressure  on  AB  can  now  be  found 
and  a little  rearrangement  gives  the  value 


P = c cot  9 


| tan2  ( f + | ) 


sln(pQ-  p) 
sln(P0+H) 


(20) 


When  9=0  (Prandtl-Reuss  material),  p Is  also  zero  and  (20) 
takes  the  form 


P = 2c(l  - cos  0O), 


agreeing  with  the  value  obtained  by  Prager  [8]#  For  small  values 
of  0q,  expression  (20)  is  approximately 

2 

P = c0Q  cos  9 , 


and  it  follows  that  for  small  values  of  the  ratio  P/c 

decreases  as  the  angle  of  friction  9 is  increased.  Fig.  4a 
shows  the  variation  of  P/c  with  the  angle  of  the  wedge  for 
9 = 0°,  20°  and  40°.  The  limiting  case  9=  90°  is  also 

shewn  in  the  diagram* 

When  = n/2,  the  expressions  (17)  and  (20)  have  the 
common  value  2c  tan  (itA  ♦ 9/2),  and  the  discontinuous  solution 
is  a continuation  of  the  continuous  solution*  The  full  lines  in 
Fig.  5b  shows  the  variation  in  P/c  with  the  wedge  angle  0O 
as  furnished  by  the  two  solutions  for  angles  of  friction  of  0°, 

20  and  40°*  The  discontinuous  solution  also  satisfies  all  the 
stress  conditions  when  the  angle  of  the  wedge  is  greater  than  n/2, 
and  the  stress  field,  Fig.  3c,  is  an  alternative  solution  to  the 
continuous  solution  of  Fig.  3a.  However,  it  will  be  seen  that  a 
velocity  field  cannot  be  associated  with  the  stress  field  of  Fig. 
3c,  so  that  the  discontinuous  solution  is  not  acceptable  physically 
for  wedge  angles  greater  than  a right  angle.  The  pressure  P as 


-fa****'—- 

4 ys 
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finished  by  the  discontinuous  solution  for  obtuse  angled  wedges  is 
again  given  by  (20)  and  this  value  is  lower  than  the  value  (17) 
obtained  from  the  continuous  solution.  The  broken  lines  in  Fig. 

4b  give  the  values  of  P/c  for  varying  £*/2)  for  angles 

of  friction  of  0°,  20°  and  40°. 

It  should  perhaps  be  pointed  out  that  the  line  of 
discontinuity  BC  in  Fig.  3c  is  not  the  line  joining  th9  point 
B in  Fig.  3a  to  the  point  of  intersection  of  the  lines  AC,  ED 
except  when  the  angle  of  friction  is  zero. 

4.  Wedge  under  unilateral  pressure.  Velocity  solutions. 

Without  restricting  the  generality  of  the  discussion, 
the  boundary  condition  on  the  velocity  field  in  the  problem  of  the 
previous  section  may  be  taken  to  be  that  the  normal  velocity  must 
have  a given  distribution  on  the  loaded  part  AB  of  the  wedge. 

The  velocity  in  the  plastic  region  has  to  be  determined  from  this 
boundary  condition  and  the  boundary  condition  at  the  plastic  rigid 
boundary.  The  given  normal  velocity  along  AB  must  be  such  that 
the  resulting  velocity  field  has  non-negative  dilatation  every- 
where and  also,  in  the  discontinuous  stress  solutions,  the  lines  of 
stress  discontinuity  must  behave  as  inextensible  filaments* 

Considering  first  the  continuous  stress  field  of  Fig;  3a, 
the  line  ACDE,  which  is  a first  failure  line,  separates  the  region 
of  plastic  flow  from  the  material  which  remains  at  rest.  The 
velocity  along  ACDE  must,  therefore,  be  inclined  at  an  angle 
9 to  ACDE  so  that  v5  is  aero  along  ACDE.  Since  the  second 
failure  lines  ar«  straight  throughout  the  plastic  region  it  follows 
that  Vg  is  zero  everywhere.  The  value  of  v^  along  AB  is 
known  from  the  giver*  normal  velocity  of  AB  and  the  fact  that 
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v2  Is  zero  along  AB.  The  variation  of  v1  along  the  first 
failure  lln*s  is  given  by  the  first  of  equations  (10)  and  knowing 
v along  AB  it  is  a simple  matter  to  obtain  throughout 

«L 

the  plastic  region.  The  condition  of  non-negative  dilatation 
will  be  satisfied  provided  that  the  given  normal  velocity  of  a 
point  on  AB  is  a non-decreasing  function  of  the  distance  of  the 
point  from  A. 

Considering  now  the  discontinuous  stress  field  for  the 
acute  angled  wedge , Fig.  3b,  the  velocity  field  must  be  such  that 
the  discontinuity  line  BC  moves  as  an  inextenslble  filament. 

Since  the  material  below  the  plastic  rigid  boundary  ACD  remains 
at  rest,  this  implies  that  the  velocity  component  along  BC  is 
zero  for  points  on  BC.  In  determining  the  velocity  field  it  is 
convenient  to  use  oblique  axes  a,0  having  A as  origin  and 
directed  along  the  first  and  second  failure  lines  at  A,  Fig.  5* 
The  normal  velocity  f(q)  is  prescribed  along  AB,  the  line  A3 
being  given  by  a * p = tj. 

The  velocity  field  in  tho  region  ABC  is  determined  in  a 
very  similar  manner  to  that  used  by  Lee  [93  for  the  same  problem 
in  a Prandtl-Reuas  material.  It  is  unnecessary,  therefore,  to 
give  the  details.  The  value  of  v^  in  regions  £ and  £ of 
Fig.  5 is  found  to  be 

vx(p)  = 2 cos(  \ - | ) f(p). 

The  velocity  projection  Vg  is  zero  in  region  £ because  of  the 
boundary  condition  along  AC,  while  in  regions  £ and  £ it  is 
given  by 

v2(a)  = -2  cos(  £ - | ) tf(  ^ ), 
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where  t » sin  sec(\|>-9),  \J>  is  the  angle  BCF  and  where  the 
length  of  AC  is  taken  to  be  unity.  For  other  points  in  the 
region  ABC,  and  v2  are  given  in  the  form  of  finite  series. 

As  the  apex  B is  approached,  the  limiting  exnressions  take  the 
form 


vx(p)  = 2cos<§  - f)  f(p)-t[2cos(£  - |)  f(£§i)  - 

-t{2cos(£  - |)  f(  ^ «—}], 

v2(a)  «=  -t[2cos(^  - f(2^i)  -t{2cos(^  - f^-^.-. j~ — J] 


The  continued  fractions  which  are  the  arguments  of  the  function 
f in  these  expressions  represent  the  coordinates  of  points  on 
AB  after  successive  reflections  at  BC  and  AB  along  lines 
parallel  to  the  oblique  axes  at  A. 

It  can  readily  be  shown  that  the  velocities  v^,  v2 
approach  limits  at  the  vertex  B.  Also  the  limits  are  such  that 
«he  velocity  of  the  plastic  field  near  B approaches  the  velocity 
of  the  non-plastic  filament  BC  at  B,  this  velocity  being  unique- 
ly determined  by  the  velocity  condition  on  BC  and  the  given 
normal  velocity  of  AB  at  B, 

The  velocities  of  the  points  on  the  filament  BC  are 
now  completely  determined.  The  velocity  of  a point  P in  the 
region  BCD  is  determined  from  the  velocity  of  the  filament  BC 
and  the  boundary  condition  at  the  plastic  rigid  boundary  CD,  The 
velocity  at  P depends  only  upon  the  velocities  of  the  points  X, 

Y where  the  failure  lines  through  P meet  the  line  BC  (see  Fig. 
5).  is  constant  along  PX,  v2  is  constant  along  PY  and 
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and  Vg  are  knovn  at  the  points  X,Y.  If  ? lies  in  the 
region  CDE  therr  v2  is  zero. 

The  function  f ( n)  which  specifies  tho  distribution  of 
normal  velocity  along  AB  can  not  be  chosen  in  a completely 
arbitrary  manner.  Certain  restrictions  must  be  placed  upon 
f(r))  in  order  that  the  dilatation  be  positive  everywhere  in  tho 
plastic  region.  The  restrictions  are  most  easily  obtained  in 
terras  of  the  velocity  X(^)  of  the  non-plastic  filament  at  the 
point  distance  t from  C.  is  continuous  since  the  filament 

is  non-plastic  and  also  x(0)  “0  from  the  condition  at  C. 
Omitting  tho  details  of  the  analysis,  it  is  found  that  the  dilata- 
tion will  be  positive  everywhere  provided  that  x*^)  2.  0 1111(1 

X'Ot)  is  a non-decreasing  function.  Those  conditions  imply  that 
the  rate  of  change  of  curvature  of  tho  filament  increases  towards 
the  vertox.  Thus  tho  wedge  behaves  in  a similar  fashion  to  a 
beam  which  is  bent  by  uniform  loading. 

The  restrictions  which  havo  to  be  placed  upon  f (n ) are 
obtained  by  exprossing  f(q)  in  terms  of  tho  restricted  function 

xa>. 

In  tho  discontinuous  stress  solution  of  Fig.  3°,  the 
velocity  conditions  for  the  problem  of  the  acute  angled  wedge 
apply  in  this  case  also.  Tho  plastic  rigid  boundary  AC  enforces 
v2  = 0 along  AC  and  it  follows  that  * 0 throughout  region 
ABC  of  Fig.  3c,  since  the  failure  lines  are  straight.  The  condi- 
tion that  the  velooity  component  along  BC  is  zero  for  points  on 
BC  then  requires  v^  to  be  zero  on  BC  and  in  consequonce  v-^ 
must  be  zero  throughout  tho  region  ABC.  Thus  the  velocity  field 
obtained  roquiros  f(n)  " 0 along  AB  and  the  solution  does  not 


All-81 


16 


constitute  plastic  flow.  The  discontinuous  stress  solution  for 


the  obtuse  angled  wedge  is  therefore  physically  inadmissible. 


A further  example  of  discontinuous  stress  fields  is 


shown  in  Fig.  6a.  The  trapezoid  ABCD  is  in  a plastic  state  of 


stress  due  to  the  normal  pressures  P,Q  on  the  parallel  sides 
AB,CD,  the  sides  AD,BC  being  free  from  applied  traction.  The 
lines  AO;  BO,  CO,  DO  are  lines  of  stress  discontinuity  separating 
the  regions  a,  £,  c,  d of  constant  stress,  where  the  angles  marked 
Y,  b in  the  figure  remain  to  be  determined. 


The  stress  boundary  conditions  require  that 


®a  " *(!  + I5' 


Ojj  «=  -(jj  + !£)  + X. 

®c  “ "C5  + f5*' 


where  2X  is  the  angle  between  the  sides  AD,  BC  and  where  the 


subscripts  refer  to  the  regions  marked  a,  b,  c in  the  figure, 
With  these  values  of  0,  the  jump  condition  (13)  at  the  dis- 


continuity lines  BO,  CO  show  that 


cos  (2Y  + X)  ° sin  9 cos  X, 
cos  (2  b - X)  **  sin  9 003  X, 


and  the  relevant  roots  of  these  equations  are 


v v X 

r-  5 - 5, 


& » £ + * , 
2 2 1 


cos  v m sin  9 cos  x> 


Oi*i§  . 
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The  values  of  p In  the  regions  3 and  c can  be 
deduced  from  the  known  value  of  p in  region  £ and  the  Jump 
conditions  (12)  at  the  lines  BO,  CO.  The  pressures  P,  Q can 
then  be  obtained  and  it  is  found  that 

P = 0 cot  t{  -.an2<tf  + f)  - x} 

Q = 0 cot  v { tan2  + f ) - 1} 

As  the  angle  of  friction  9 tends  to  zero  these  expressions  tend 
to  the  values 

PQ  « 2c(l+sin  X),  Q0  * 2c(l-sin  x)> 

agreeing  with  the  values  obtained  by  Winzer  and  Carrier  [ 12]  for 
a Prandtl-Reuss  material. 

It  can  be  shown  that  a velocity  field  can  not  be 
associated  with  the  stress  field  of  Fig.  6a.  This  follows  in  the 
same  way  in  which  it  was  shown  that  a velocity  field  could  not 
be  associated  with  the  stress  field  of  Fig.  3o.  The  stress  field 
is  therefore  physically  inadmissible. 

If  the  sides  AD,  BC  of  the  trapezoid  are  continued 
beyond  the  points  D,  C a truncated  wedge  is  obtained.  A solu- 
tion corresponding  to  a uniform  pressure  on  the  top  surface  of 
the  wedge  has  been  given  by  Prandtl  [9]  and  the  stress  field  is 
outlined  in  Fig.  6b,  The  value  of  the  normal  pressure  P*  is 
in  this  caso 

P»  • c oot  9 ftan2(£  + 3f)  exp(2x  tan  9)  - l]  . 

and  this  value  is  greater  than  the  value  P given  by  equation 
(22).  A velocity  flold  oan  be  associated  with  the  stress  field 


(22) 
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of  Fig.  6b  provided  that  the  given  normal  velooity  of  the  points 
on  AB  is  constant  or  a monotonic  function  of  the  distance  from 
A. 

6.  Integration  of  the. jplan^train  equations . 

In  obtaining  solutions  to  plasticity  problems  in  soil 
mechanics  the  stress  distribution  can  not  always  be  built  up 
from  the  simple  failure  patterns  of  the  Rankine  zones  and  the 
zone  of  radial  shear.  In  complox  problems  involving  fields  of 
varying  stress  it  is  usual  to  obtain  tho  stresses  by  approximate 
numerical  methods  [3].  In  this  section  exact  expressions  will  be 
dorived  for  tho  coordinates  and  curvatures  of  the  failure  lines 
and  tho  velocity  components  at  any  point  of  a stress  field  in  terms 
of  the  boundary  valuos.  However,  the  expressions  involve  integrals 
which  can  be  evaluated  in  finite  form  only  whon  the  boundairy 
conditions  are  particularly  simple.  For  this  reason  tho  analysis 
is  of  limitod  importance  and,  as  in  tho  theory  of  a Prandtl-Rouss 
material,  will  mainly  be  of  value  in  checking  tho  accuracy  of  the 
usual  approximate  methods. 

In  Fig,  7,  the  x,  y-axes  are  taken  at  a point  0 along 
and  normal  to  the  direction  of  tho  first  failure  line  at  0.  The 
first  and  second  failure  linos  0A  and  OB  will  be  called  the 
a base-line  and  the  p base-line  respectively.  The  curvilinear 
coordinates  of  a point  P are  defined  by  the  pair  of  quantities 
(a,p);  where  a is  the  value  of  0 at  the  point  A where  the 
second  failure  line  through  P meets  the  a base-line  0A,  and 
where  0 is  the  value  of  0 at  the  point  B where  the  first 
failure  line  through  P meets  the  0 base-line  CB.  The  corre- 
spondence between  (a,P)  and  the  coordinates  (x,y)  of  P is 
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one  to  one  except  when  one  of  the  families  of  failure  lines  is 
straight.  It  follows  from  the  definition  that  a and  0 are 
constant  along  the  second  and  first  failure  lines  respectively, 
and  the  base-lines  are  given  by  0=0  and  u = 0. 

With  equations  (5)  it  is  not  difficult  to  show  that 


0 = a 0 , 


1 

2 


cot  9 log  p = 0 


a 


1 

2 


cot  9 log  P0, 


(22) 


where  0 and  p are  evaluated  at  the  point  P and  p^  is 
the  value  of  p at  the  origin  0.  The  first  of  equations  (23) 
is  a statement  of  what  is  usually  known  as  Hencky's  first  theorem 
[13],  which  is  expressed  analytically  by 


6P  - «B  = «A  - V (2k) 

The  theorem  was  first  stated  by  Massau  [4]  for  the  case  of  a 
cohesionless  soil.  Equation  (23)  can  be  re-arranged  to  give 


tions 


(25) 


20  = - 2 cot  9 pQ  + ( \ cot  9 log  p + 0).  J 
Coordinates  x,y  are  introduced  by  means  of  the  equa- 


x = x cos(©+9)  + y sin(©+<p),  y = -x  sin  © + y cos  0, 
x = (x  cos  © - y sin(©+<p)}  sec  9,  y = [x  sin©  + y cos(©+9>}  sec  9. 

* 

(26) 

The  geometrical  significance  of  the  coordinates  (x,y)  is  shown 
in  Fig.  8;  x and  y are  the  distances  of  the  origin  0 from 
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the  tangents  to  the  second  and  first  failure  lines  at  P.  The 
coordinates  x,y  have  been  used  in  the  theory  of  the  Prandtl- 
Reuss  material  [14],  On  a first  failure  line  (on  which  dy  = 
tan  9 dx)  equations  (26)  give 

dy  = - sin  9 dx  + cos  8 dy  - (x  cos  8 + y sin  8)  d8 
= - (x  sec  9 - y tan  9)  d8, 

and  therefore 

dy  +(x  sec  9 -y  tan  9M8  = 0 along  a first  failure  lino.  (27) 
Similarly, 

dx  -(y  sec  9 -x  tan  9M8  = 0 along  a second  failure  line.  (28) 

The  radii  of  curvature  R,S  of  the  first  and  second 
failure  lines  »t  the  point  P are  defined  by 


(29) 


where  s^  and  Sg  are  the  arc-lengths  along  the  first  and  second 
failure  lines.  Pig.  9 shows  the  intersection  of  two  neighboring 
first  failure  lines  with  two  neighboring  second  failure  lines. 

Tho  small  angles  8g  - 8^,  8^  - 8^  will  bo  denoted  by  A8^,A82 

respectively.  In  ; he  figure,  A8^  is  positive  while  A82  is 
negative.  It  follows  from  the  definitions  (29)  that 


RA81  =Asx,  SA82  = -As2 


(30) 


and  therefore 


0 

0 s„ 


(3D 


! 


All-Sl 

From  the  geometry  of  Fig.  9 and  equations  (30), 


21 


DC  = As^  - A9^  sec  9 A s^  - A©2  tan  9 A s., 

= As1  - A 9^82 (sec  9 - ■§  tan  9), 

where  the  sign  of  AOg  has  been  taken  into  account.  Hence 

(As  ) = -A9,  (sec  9 - S tan  9).  Equation  (24)  shows  that 
os  1 1 S 

A 9-.  is  a constant  and  it  can  be  taken  outside  the  operator  on 
the  left  hand  side  of  equation  (31),  and  we  have  finally 


Similarly 


AS  = -(sec  9 - •§  tan  9). 
8 s2  Y S 


= -(sec  9 - | tan  9) . J 


(32) 


These  equations  give  the  variation  of  the  radii  of 
curvature  along  the  failure  lines.  The  corresponding  equations 
for  a cohesionless  heavy  soil  have  been  deduced  previously  by 
Massau[4].  When  9=0  equations  (32)  become 


AS  = 
9S2 


-1 


» 


= -1. 


and  this  is  an  analytical  statement  of  Hencky’s  second  theorem 

for  a Prandtl-Reuss  material. 

The  substitution  of  RA9,  and  -SAO^  for  the  arc- 

l 2 

lengths  as^  and  As2  gives  an  alternative  form  of  equations 
(32), 


dS  +(R  sec  9 - S tan  9M9  = 0 
dR  -( S sec  9 - R tan  9M9  = 0 


along  a first  failure  line, 
along  a second  failure  line. 


(33) 


1 


i 


K 
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On  a first  failure  line  dO  = da  and  dO  = p cn  a 
second  failure  Ifne  so  that  equations  (10)  can  be  written 


An  alternative  form  of  these  equations  is 


Satv  1exp{(p-a)  tan  9}]-  secqi  v2  exp^(p-a)  tancpjc 

g£-[v2exp{(p-a)  tan  9}]+  sec9V^  exp{(p-a)  tan9}  * 

In  the  same  way,  equations  (27)  and  (28)  can  be  replaced 
JL  [ y exp  {(jp -a)  tan  9]]+  sec  9 xexpj(p-a)  tan  9}=  0, 

^[x  exp(<j>-o)  tan  9}]  - see  9 y exp^O-a)  tan  9}=  0, 
and  equations  (33)  become 

•&[S  exp {(p  -a)  tan  9]  3+  sec  9 R exp  {(p -a  Kan  9}  = 0,  ] 

^£r  exp{(p-a)tan  9}]-  sec  9 S exp{(p-a)tan  9}  = 0.  J 

With  equations  (23),  the  exponential  term  in  equations  (34-),  (35), 
(36)  can  be  written 

expf(p -a)tan  9}  = exp  j log -E-)  = / 

12  V IV' 

It  can  easily  be  deduced  from  equations  (34),  (35),  (36) 
that  each  of  the  quantities  v^,  v2,  x,  y,  R,  S multiplied  by 
exp{Q)-a)  tan9^  satisfies  the  differential  equation 


by 

) 

V (35) 


(36) 
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2 

-Hz  + sec2  9 f = 0.  (37) 

The  method  of  solution  of  this  equation,  which  is  one  form  of  the 
equation  of  telegraphy,  depends  upon  the  manner  in  which  the 
boundary  conditions  on  the  function  f are  prescribed.  The 
analytic  solution  for  one  type  of  boundary  condition  is  given 
below.  The  solutions  for  other  types  follow  closely  the  analytic 
solutions  for  the  corresponding  problems  in  a Prandtl-Reuss 
material  [10  ]« 

We  suppose  that  the  stress  components  are  given  along 
a curve  T which  is  not  a failure  line.  If  only  the  normal 
and  tangential  stresses  are  given  on  X,  two  values  for  the 
interior  stress  component  can  be  obtained  from  the  yield  condition 
and  the  correct  value  must  be  determined  from  the  conditions  of 
the  problem  in  question.  The  stressos  on  X define  the  variation 
cf  p and  0 along  X,  and,  from  the  theory  of  characteristics, 
these  boundary  values  define  the  stress  field  in  the  region  DEC 
enclosed  by  X and  the  two  failure  lines,  one  from  each  family, 
passing  through  the  end  points  D,E  of  X and  meeting  at  a 
point  C.  In  Fig.  10a,  the  failure  line  DC  is  taken  to  belong 
to  the  first  family  of  failure  lines.  The  values  of  a and  0 
along  X can  be  obtained  from  equations  (25)  and  Fig.  10b  shows 
the  configuration  in  the  a{5  plane. 

The  values  of  R and  S along  X are  obtained  as 


follows.  If  s is  the  arc-length  along  X from  the  point  D 
then 


(38) 


where  is  the  angle  between  the  first  failure  direction  and 
the  tangent  to  Z,  Hence,  from  equation  (29), 


-22  = cos(\|>-9)  sec  <p/R  - sin  \|>  sec  9/S  , 


and  also 


JL[JL  cot  9 log  p ] = - cos(>l>-9)  sec  9 /R  - sin  tl>  sec  9 /S 
0 s 2 

fr<-n  equations  (>)  and  ( 29) • R and  S are  therefore  given  at 
points  on  Z by  the  relations 

(39) 


The  derivatives  of  R and  S with  respect  to  p and  a respective- 

ly are  known  along  2 from  equations  (36)  when  R and  S are 
known  along  Z « The  other  derivatives  can  then  be  calculated  from 
the  equation 

££  = {costy-9 ) sec  9/h]  - (sin  sec  9/S  } §£  , (l+o) 

which  follows  from  equation  (38)  and  the  definitions  of  R,  S. 

The  boundary  values  of  (x,  y)  are  given  immediately 
by  equations  (26)  so  that  the  boundary  values  of  the  quantities 
x exp  £(p  -o)  tan  9],  y sxp^(P-a)  tan  9}  can  be  found.  The 
derivatives  of  these  quantities  at  points  on  Z are  calculable 
from  equations  (35)  and  (40). 

We  suppose  also  that  v^  and  v2  are  specified  on  the 
curve  2 • The  boundary  values  assigned  to  v^,  Vg  must  be 
compatible  with  the  stress  boundary  conditions  on  Z in  order 


1 

R 


1 

s 


-!-££££—  -2- [I  cot 9 
2 cosfy  -9  ) ® 3 


log  p - 0], 


= - -k c i oot  f log  p + 9] 
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to  ensure  that  the  resulting  velocity  solution  has  non-negative 
dilatation  everywhere.  The  quantities  v^exp^p-cOtan  9}  , 
v2exp((p-a)  tan  9}  are  then  known  on  Z and  the  derivatives  of 
these  quantities  at  points  on  S can  be  found  from  equations  (34) 
and  (40).  This  is  equivalent  to  specifying  a member  of  the  pair 
(v  ,v2)  and  one  of  its  derivatives  since  the  other  member  of  the 
pair  can  be  determined  from  equations  (34)  and  (4o). 

Thus  the  problem  of  determining  the  quantities  (v1>v2)> 
(x,y)  and  (R,S)  in  the  region  DEC  is  reduced  to  the  problem 
of  determining  a function  f(a,  {))  which  satisfies  equation  (37) 
and  which  has  prescribed  values  for  f and  one  of  its  derivatives 
on  the  curve  T , The  other  derivatives  at  points  on  2)  can  be 
obtained  from  equation  (40).  Let  P(a,b)  be  a point  in  the  region 
DEC  and  lot  AP  and  BP  be  the  failure  lines  through  P meeting 
the  curve  Z at  A and  B.  Employing  the  method  of  Riemann,  a 
particular  solution  of  equation  (37)  with  certain  properties  to  be 
defined  later  will  be  denoted  by  G(a,(J).  It  follows  that  the 
expression 


is  a perfect 

! 

! 

1 

1 

I 

t 

The  integral 
is  therefore 
curve  APBA 


(4l) 


differential  since 
G 


— = -Of  sec2  <p  = t&L  . 


Oa90  T ‘ 8aep 

of  the  expression  (4l)  taken  around  any  closod  curve 
zero  and  if  the  expression  is  integrated  around  tho 
there  results 


I 

1 

! 


— T 
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rP 


0 = 


fB 

(f  ap-0#^  4P  + ! (G  f & *» 

^ A 


j, 


6 a Oa 


(42) 


[(Gsf- f 88>  * + (f  $-  Gt^>  w) 


B 


In  order  tc  evaluate  the  first  two  integrals  explicitly,  the 
function  G(a,  £)  is  now  defined  to  be  such  that  G = 1 on  AP 
and  BP.  This  implies  that  = 0 on  AP  and  = 0 on  BP 
and,  after  some  re-arrangement,  equation  (39)  gives 


vi<w  * i Jb(<g  fa  -f  <f  $ - G ■§£)  • 0*> 

It  can  be  shown  that  the  particular  function  G is  given  by 
G(c,p)  = J0[2  sec  9{(a-a)(b-p)}*], 


where  Jq(^)  is  the  Bessel  function  of  order  zero, 

j * 

Jo  + -f  * Jo  = 0.  J0<°>  * 1.  Jo<°)  = 0. 

Equation  (*+3)  therefore  expresses  the  value  of  f at  the  point  P 
in  terms  of  f and  its  derivatives  along  A3,  f can  be  replaced 
by  any  one  of  the  quantities  v^,  v2 , x,  y,  R,  S multiplied  by 
exp{(0-u)  tan  9}, 
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FIG.  4a  PRESSURE  ON  FACE  OF  ACUTE  ANGLED  WEDGE 
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FIG.  4b.  PRESSURE  ON  WEDGE  FACE  AS  GIVEN  BY  THE  DISCONTINUOUS 

AND  CONTINUOUS  SOLUTIONS 
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RG.  ta  FAILURE  LINES  AND  BOUNDARY  CURVE  X IN  BOUNDARY  VAUUf 

PROBLEM  CONSIDERED 


